Energy 

Conversion & 
Management 

Energy Conversion & Management 41 (2000) 235-248 - 

www.elsevier.com/locate/enconman 

A comparative performance analysis of irreversible Carnot 
heat engines under maximum power density and maximum 

power conditions 

Ali Kodal a , Bahri Sahin b ’*, Tamer Yilmaz b 

a Department of Aeronautical Engineering, Istanbul Technical University, Maslak, 80626, Istanbul, Turkey 
b Department of Naval Architecture, Yildiz Technical University, Besiktas, 80750, Istanbul, Turkey 

Received 3 November 1998; accepted 3 June 1999 



PERGAMON 


Abstract 

This paper reports the finite time thermodynamic optimization based on the maximum power density 
criterion for an irreversible Carnot heat engine model which includes three types of irreversibilities: finite 
rate heat transfer, heat leakage and internal irreversibility. The obtained results are compared with those 
results obtained by using the maximum power criterion. The design parameters under the optimal 
conditions have been derived analytically and the effects of the irreversibilities on the general and 
optimal performances are investigated. The results showed that the design parameters at maximum 
power density lead to smaller and more efficient heat engines. It is also seen that the irreversibilities have 
a greater influence on the performances at maximum power density conditions with respect to the ones 
at maximum power conditions. Also in this study, the optimal conductance allocation parameter is 
investigated at both maximum power density and maximum power conditions by assuming a 
constrained total thermal conductance in the case when there is no heat leakage. The relation between 
the conductance allocation parameter and the thermal efficiencies at maximum power density and 
maximum power is also investigated. The obtained results generalize the result of previous studies on 
this subject and provide guidance to optimal design in terms of power, thermal efficiency and engine 
sizes for real heat engines. © 1999 Elsevier Science Ltd. All rights reserved. 
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Nomenclature 

k specific heat ratio 

m mass of working fluid 

m mass flow rate 

P pressure 

Q rate of heat transfer 

Q m total heat rate transferred from hot reservoir 

2lt total heat rate transferred to cold reservoir 

R ideal gas constant 

R/ss internal irreversibility parameter 

S entropy 

T temperature 

V 4 maximum volume in cycle 

W power generated from heat engine 

.v conductance allocation parameter (a/2) 

ot,P,y thermal conductance of heat source, heat sink and heat leakage sides, respectively 
t] thermal efficiency 

X overall thermal conductance when no heat leakage (a + /J) 

r ratio of heat-sink temperature to heat-source temperature 

Subscripts 
C Carnot 

d density 

H heat source 

L heat sink 

LK heat leakage 

max maximum 

min minimum 

mp maximum power conditions 

opt optimum 

X warm working fluid 

Y cold working fluid 

Superscripts 

* maximum power density conditions 


1. Introduction 

Power optimization studies of heat engines using finite time thermodynamics analysis were 
started by Chambadal [1] and Novikov [2] in 1957 and 1958 and were continued by Curzon 
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and Ahlborn [3] in 1975. They extended the reversible Carnot cycle to the endoreversible cycle 
by taking the irreversibility of finite-time heat transfer into account. By introducing finite time 
processes, a new field: ‘finite time thermodynamics’ was born, and thus, more realistic limits 
have been adapted for real heat engine performances. During the last decade, many 
optimization studies for heat engines based on endoreversible and irreversible models have 
been performed by considering finite time and finite size constraints under various heat transfer 
modes, mainly linear and non-linear ones. The reader who is interested in these optimization 
works could refer to a literature survey written by Bejan [4], Usually, in these studies, the 
power output, thermal efficiency, entropy generation and the ecological benefit were chosen for 
optimization criteria. However, the performance analyses based on the above optimization 
criteria do not take the effect of engine sizes related to the investment cost into account. In 
order to include the effects of engine size in the performance analysis, Sahin et al. [5] 
introduced a new optimization criterion called the maximum power density (MPD) analysis. 
Using the MPD criterion, they investigated optimal performance conditions for reversible [5] 
and irreversible [6] non-regenerative Joule-Brayton heat engines. In their study, they 
maximized the power density (the ratio of power to the maximum specific volume in the cycle) 
and found design parameters at MPD conditions, which lead to smaller and more efficient 
Joule-Brayton engines than those engines working at maximum power (MP) conditions. Erbay 
et al. [7], Erbay and Yavuz [8], Chen et al. [9] and Medina et al. [10] applied the MPD 
technique to the Ericsson, Stirling, Atkinson and regenerative Joule-Brayton engines, 
respectively. In their analyses, they discussed the advantages of the MPD performance 
conditions in comparison to the MP conditions in terms of thermal efficiency and engine sizes. 
Sahin et al. [11] applied the MPD technique to the endoreversible Carnot heat engine which 
can be considered as a theoretical comparison standard for all real heat engines in finite time 
thermodynamics and thus generalized the endoreversible MPD analyses results. They also 
demonstrated that the thermal efficiency at MPD conditions varies between 1 — Jr and [1 — 
2t/(1 + t)] for 0<v< 1 where t is the ratio of heat sink to heat source temperature and x is the 
conductance allocation parameter. 

In this paper, the MPD performance analysis performed by Sahin et al. [11] for the 
endoreversible Carnot heat engine model is extended to an irreversible Carnot heat engine 
model by considering internal irreversibility and heat leakage effects. In this context, the 
optimal performance and design parameters under MPD conditions are investigated. The 
obtained results are comparatively discussed with respect to the results obtained by using the 
MP performance criterion. 


2. The theoretical model 

The considered irreversible Carnot heat engine model which includes heat leakage, finite time 
heat transfer and internal irreversibilities and its T-S diagram are shown in Fig. 1. The rate of 
heat leakage (2 lk from the hot reservoir at temperature r H to the cold reservoir at temperature 
7) with thermal conductance y is given by 


£?lk = 7<Th - 71). 


( 1 ) 
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Fig. 1. Irreversible Carnot heat engine model and its T-S diagram. 


The heat rate Q H from the hot reservoir at temperature T H to the heat engine with thermal 
conductance a and the heat rate Q L from the heat engine to the cold reservoir at temperature 


T\ with thermal conductance [>. are respectively, 

Qh = a(7 h - T X ), (2) 

Ql = p(T Y - T l ). (3) 

Then, the total heat rate g HT transferred from the hot reservoir is 

Qm — Qh + 2lk = a (2"H — T x ) + y (T H — TL) (4) 

and the total heat rate Q LT transferred to the cold reservoir is 

Qlt = Ql + 2lk - P(Ty ~ TL) + y(T H - TL), (5) 

where T x and T Y are the warm and cold working fluid temperatures, respectively. The power 
produced by the engine according to the first law is 

W = Qm - Qlt = Qh~Ql = «(T h - T x ) - p(T Y - TL). (6) 

The power density defined as the ratio of power to the maximum volume in the cycle then 
takes the form 

Wd = [a(T H - T x ) - p(T Y - TL)] / V 4 . (7) 

Formally, the second law for an irreversible cycle requires that 
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SQ = Qh_Ql <0 
T T x T y 


One can rewrite the inequality in Eq. (8) as 
^ - Ras^J = 0 with 0 < r as < 1 

tx 1 Y 

where the internal irreversibility parameter Ras is defined as 


Ras — 


S3 ~ S 2 

s 4 -s 1 


which characterizes the degree of internal irreversibility. By substituting Eqs. (2) and (3) in Eq. 
(9), we have 

a(E h - T x )/T x - RasP(T y - T L )/T Y = 0. (11) 

Assuming an ideal gas, the maximum volume in the cycle V 4 can be written as 

V 4 = mRT Y /P min , (12) 

where m is the mass of the working fluid and R is the ideal gas constant. In the analysis, the 
minimum pressure (P Trnn ) in the cycle is taken to be constant. The power density then becomes 

Wd = (P mm /niR)[^(T H - T x ) - P(T y - T L )]/T Y . (13) 

One can maximize the power density given in Eq. (13) with respect to T x and T Y by using Eq. 
(11). The results are 

T*x = W(a + M/(* + PRas), (14) 


T y — PRasTJ a + PRas ~ <Xy/(a + PR Y s)/{u. + ~jh) , 


where z = T L /T H . Substituting Eqs. (14) and (15) into Eq. (13), the MPD can be found as 

r _ _-i2 

= (y.P mm /mR) y/a + PR& s - Jet, + pz /PzR as . 

The thermal efficiency, 

— 1 ~ Q lt /Q ht 

can be written, after some manipulations as 

0 = [1 - Ty/CRasEx)]/] 1 + [y(l - t)(« + PRas)T y /T x ]/[^PRas(T y /T x - r)] j. 


Using Eqs. (14) and (15) in Eq. (18), the efficiency at MPD becomes 
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>f = o/i 1 + 


■ y(l - t)/ a(l - 7( a + jSx)/(a + [IRas)) 


(19) 


where ;y* =0 is the efficiency at MPD when there is no heat leakage, and it can be obtained as 


n * y= o = 1 - j6t/(V (a + [h)(y. + fiR&s) - a)- 


( 20 ) 


Similarly, one can maximize the power given in Eq. (6) with respect to T x and T Y using the 
condition given in Eq. (11) and obtain 


T Xmp = (« 7 h + IWXasT l T h )/(v + pR AS ) 

T Ymp = (pRT L + ol-Jr as2l7h)/(« + PRas)- 
Substituting Eqs. (21) and (22) into Eq. (6), the maximum power can be found as 

= y-fiT H (/R^s - \/r) /(a + fiRtss) 


W — 

rv max — 

and the thermal efficiency at maximum power is 

i - V^JRas 


n mp = 


1 + [y(l - t)(0C + PRas)]/[uPRas(1 - \/ T /^As)] 


( 21 ) 

( 22 ) 

(23) 

(24) 


3. Results and discussion 

The variations of the normalized power density and the normalized power with respect to 
thermal efficiency are shown for the zero heat leakage case (y = 0) in Fig. 2(a) and (b), and for 
the heat leakage case (y > 0) in Fig. 3(a) and (b), respectively. The following comments can be 
summarized from the observation of these figures: (i) for the chosen values of the parameters, 
although the maximum thermal efficiencies 0/ max ) are the same, the thermal efficiency at MPD 
(if) is greater than the thermal efficiency at MP conditions (»/ mp ), i.e, ;y max > if > )/ mp ; (ii) an 
increase in t (the ratio of heat sink to heat source temperatures) and the effect of heat leakage 
cause the thermal efficiencies »/ max , if and j/ mp to get closer to each other, that is the thermal 
efficiency advantage of the MPD conditions with respect to the MP conditions (A;/ = if — j/ mp ) 
decreases; (iii) the heat leakage effect causes the loop shaped performance curves and also 
causes decreased values of the optimal thermal efficiencies, e.g. for r —0.1 and R AY = 0.8 
values when there is no heat leakage (y = 0), ;/ max — 0.88, if = 0.75 and t] — 0.64 and when 
there is heat leakage (y/a = 0.05), ;/ max = 0.62, if = 0.60 and ;/ mp — 0.55. The above 
evaluations, drawn from observation of Figs. 2 and 3 can be seen more clearly from Fig. 4(a). 
The variations of the thermal efficiencies at MPD and MP conditions (if and j/ mp ) and the 
variations of the dimensionless powers at MPD and MP conditions [W — WocT H ,tV mp — 
fffmax/aT’n] with respect to t can be seen from Fig. 4(a) and (b), respectively. In Fig. 4(a), the 
Carnot efficiency (ij c = 1 — t) is also given for comparison reasons. As can be seen from Fig. 
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X = 0.1 0.3 0.5 0.7 



*1 

(a) 


T = 0.1 0.3 0.5 0.7 



(b) 


Fig. 2. Variations of the normalized power density (a) and the normalized power (b) with respect to the thermal 
efficiency (a = [3, y = 0, Ras = 0.8). 



Fig. 3. Variations of the normalized power density (a) and the normalized power (b) with respect to the thermal 
efficiency (a = [3, y/a = 0.05, R\s = 0.8). 
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Fig. 4. Variations of the thermal efficiencies (a) and dimensionless powers (b) at MPD and MP conditions with 
respect to t, (a = [). R&s = 0.8). 


4(a) and (b), if is greater than t] mp and W* is lower than that is if the design parameters 

are selected at MPD conditions instead of MP conditions the thermal efficiency increases as 
much as A;/ — if — ;/ mp , but the power reduces by A W = W mm — W . On the other hand, as t 
increases, if and j/ mp , and W and W m ax decrease and get closer to each other, respectively and 
they become zero when x = R,\s- In order to obtain positive performance in terms of thermal 
efficiency and power for both the MPD and MP conditions, it is necessary that x < R^s- This 
case can be seen by evaluating Eqs. (16), (19), (23), (24). The variations of the dimensionless 
power density [Wj — W ( i/(aP m i n /mR)\ and the dimensionless power [W — W/(aT H )\ with 
respect to thermal efficiency are shown for various internal irreversibility parameters (R\s) in 
Fig. 5(a) and (b), and for various heat leakage values in Fig. 6(a) and (b), respectively. From 
these figures, one can observe the effects of the internal irreversibility and heat leakage on the 
global and optimal performances. As the internal irreversibility and the heat leakage increase, 
that is, as R/ss decreases and y/a increases, the global and the optimal performances gradually 
decrease and the values of the efficiencies ;/ max , if and ;/ mp get very close to each other. In Fig. 
7(a) and (b), the effects of the internal irreversibility and heat leakage on the thermal 
efficiencies at MPD and MP conditions can be seen more clearly. As can be seen from the 
figures, if is greater than j/ mp for all values of R&s and y/a. However, the advantage of if with 
respect to »/ mp (A)/ = if — j/ mp ) decreases as the internal irreversibility and heat leakage increase. 
Therefore, the reducing effects of the internal irreversibility and heat leakage on if is greater in 
comparison to j/ mp . Thus, it is more important to take precautions to reduce internal 
irreversibility and heat leakage in the design of a heat engine working at MPD conditions with 
respect to the one working at MP conditions. 

In order to determine optimal thermal conductances, which are measures of heat exchanger 
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Fig. 5. Variations of the dimensionless power density (a) and the dimensionless power (b) with respect to the 
thermal efficiency, for various internal irreversibility parameters, (r = 0.2, a = /?, y/a = 0.05, for the solid lines 
y = 0). 



Fig. 6. Variations of the dimensionless power density (a) and the dimensionless power (b) with respect to the 
thermal efficiency, for various heat leakage values, (t = 0.2, a. = ft, R& s = 0.8). 
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Fig. 7. Effects of the internal irreversibility (a) and heat leakage (b) on the thermal efficiencies at MPD and MP 
conditions, (t = 0.2, a = /?). 


sizes, we assume that the total thermal conductance is constrained, and there is no heat 
leakage, i.e 

a + ft — X = constant. (25) 

Let us define a thermal conductance allocation ratio x for the heat source side as 

a = xX, (26) 

and from Eq. (25) for the heat sink side 

/? = (!— x)X. (27) 


By using these definitions, the maximum power density given in Eq. (16) becomes, 




mR (1 — ,v)t Ras 


(28) 
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The effects of the conductance allocation parameter x on dimensionless maximum power 
density Wd mm = Wd,„„/(2P m m/mR) are shown in Fig. 8(a) and (b) for various i and R,\s values. 
As can be seen from the figures, the optimal conductance allocation parameter x which 
maximizes Eq. (28) does not much affected by x and Ras and it is approximately x* pt = 0.4. 

For the maximum power conditions, by using Eqs. (26) and (27) in Eq. (23), we have 


W 


max 


2x( 1 — x)[V Th — \/Ty\ 
x + (1 - x)R A s 


(30) 


and the thermal efficiency at MP given in Eq. (24) becomes, 
'/mp = 1 - V^S. 


(31) 


The optimal value of x which yields the maximum of the maximum power can be found by 
maximizing Eq. (30) with respect to x. The result of this optimization is 

( A ' m P)opt = n— 0.5 , i ( 32 ) 

which is independent of x and only depends on Ras- The optimal conductance allocation 
parameter at MPD is approximately constant around 0.4 (not much affected by x and Ras), 
but the optimal conductance allocation parameter at MP changes with Ras- For example, in 
the endoreversible case (Ras = 1), (x mp ) opt = 0.5 and it reduces to 0.47 for an irreversible case 
when Ras = 0.8. When we compare the optimal conductance allocation parameters at MPD 
and MP conditions, we observe that (x mp ) opt > x* pt for Ras > 0.45. The relations between the 
optimal thermal efficiencies (if and tj ) and the conductance allocation parameter can be seen 
from Eqs. (29) and (31). As is seen, t] is independent of x but if depends on x, this case can 
also be seen from Fig. 8(c). This figure and Eq. (29) indicate that, as x—>1, if approaches 1 — 
2x/(x + Ras) and for x = 0, if = i / mp = 1 - «Jx/R AS , i.e. 

1 - s/t/Ras< tf < 1 - 2x/(x + R as ) for 0<x<l. (33) 


This result generalizes the result previously shown by Sahin et al. [11]. 

The size of a heat engine can be characterized by the maximum volume in the cycle, i.e V 4 . 
For example, in reciprocating engines, the engine size is determined by the sum of the 
clearance and displacement volumes (i.e the maximum volume in the cycle). For 
turbomachinery, Sahin et al. [5] discussed that the maximum volume in the cycle characterizes 
the engine size. The ratio of the maximum volume at MPD to the one at MP can be written as 

V* (1 - x)yjR as^I-X + (1 - x)/?A.s][x + (1 - x)t] 

:. ■ . — ? " Y • ( 34 ) 

4 mp [x + ( 1 - x) ] I y/fx + d- X)/?A,][-V + (1 - x)t] - X | 
similarly the ratio of the maximum pressures at MPD and MP conditions is 
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t = 0.1 0.2 0.3 0.4 



X 

(a) 


R A s= 0.7 0.8 0.9 1.0 



(b) 



x 


(C) 

Fig. 8. Effect of the conductance allocation parameter x on dimensionless maximum power density for various r (a), 
(Ras = 0.8, y = 0); and for various Ras values (b), (r = 0.2, y = 0); and the effect of x on thermal efficiencies at 
MPD and MP (c), (r = 0.2, R AS = 0.8, y = 0). 
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p* 

1 max 


(.P max) 


mp 


= exp 


2m R 





1 

Tv 

■A mp 


\T*xTY mv 


k/l-k 


(35) 


In the derivation of Eq. (35), it is assumed that Si — S\ = S 4 — S 3 . Variations of V* 4 /(V 4 ) mp 
and i > ^ ax /(i , max ) mp are plotted with respect to t for the endoreversible case (R\s = 1) and an 
irreversible case (Ras = 0.8) in Fig. 9. As can be seen from the figure, the ratio of the 
maximum volumes in the cycle is less than unity and the ratio of the maximum pressures in the 
cycle is greater than unity. As a result, at MPD conditions engine sizes will be smaller but will 
require higher maximum pressure than the one operating at MP. On the other hand, as 
internal irreversibility increases, that is, as Ras decreases, the engine size advantage of the 
MPD conditions decreases. However, its maximum pressure disadvantage gets better. 


4. Conclusions 

A comparative performance analysis based on MPD and MP criteria has been performed for 
irreversible Carnot heat engines which include the major irreversibilities existing in real heat 
engines. The analysis showed that the thermal efficiency at MPD conditions is greater than the 
one at MP conditions, for an irreversible heat engine. However, the thermal efficiency 



Fig. 9. The variation of the ratio of the maximum volume in the cycle at MPD to that at MP and the variation of 
the ratio of maximum pressure at MPD to that at MP with respect to r for endoreversible case (Ras = b solid lines) 
and irreversible case (Ras = 0-8, dashed lines), (.v = 0.5, y = 0, T L = 300 K, 1/iiiR = 3). 
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advantage of the MPD conditions with respect to MP conditions decreases as the internal 
irreversibility and heat leakage increase. This result indicates that the reducing effect of the 
internal irreversibility and heat leakage on the thermal efficiency at MPD conditions is greater 
in comparison to the one at MP conditions. Thus, it is more important to take precautions to 
reduce internal irreversibility and heat leakage in the design of a heat engine working at MPD 
conditions with respect to the one working at MP conditions. It is found that the optimal 
conductance allocation parameter at MPD conditions (x* pt ) is approximately constant around 
0.4, and it is not much affected by z and R&s (Fig. 8(a) and (b)) but the optimal conductance- 
allocation parameter at MP conditions [(x mp ) opt ] changes with R&s (see Eq. (32)). When we 
compare the optimal conductance-allocation parameters at MPD and MP conditions, we 
observe that (x mp ) opt > x* pt for R&s > 0.45. When there is no heat leakage, the thermal 
efficiency at MP conditions is independent of the conductance-allocation parameter (see Eq. 
(31)) but the thermal efficiency at MPD conditions depends on the conductance-allocation 
parameter (Eq. (29)). Eq. (29) indicates that the thermal efficiency at MPD conditions varies 
between 1 — V z/R\s and 1 — 2 t/(t + R\s) for 0<x<l. Furthermore, it is shown that engine 
sizes designed at MPD conditions would be smaller but require higher maximum pressure than 
those operating at MP conditions. The results presented in this analysis generalize the results 
of the previous studies on this subject and may provide guidelines for determination of the 
optimal design and operating conditions of real heat engines. 
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